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Abstract. A general class of mass transport models with Q species of conserved 
mass is considered. The models are defined on a lattice with parallel discrete 
time update rules. For one-dimensional, totally asymmetric dynamics we derive 
necessary and sufficient conditions on the mass transfer dynamics under which 
the steady state factorises. We generalise the model to mass transfer on arbitrary 
lattices and present sufficient conditions for factorisation. In both cases, explicit 
results for random sequential update and continuous time limits are given. 



1. Introduction 

A wide variety of physical phenomena may be characterised by the stochastic transfer 
of one or more conserved quantities between regions of space. Examples include traffic 
flow PEl, shaken granular systems [3], sandpile dynamics @], cluster aggregation and 
fragmentation [5] and the dynamics of phase separation [5] . However, for models of 
such processes, which are devised on the basis of some putative microscopic dynamics 
in order to capture the essential elements of these systems, there exists no general 
prescription to identify the steady state, if indeed it exists at all. 

Among the models that have been proposed, particularly well studied examples 
are the zero-range process (ZRP) |H] and the asymmetric random average process 
[HI E3j- These models involve mass transfer between the sites of a lattice. Chief 
among the reasons for the widespread interest in these models is that they can be 
solved exactly in the steady state, a virtue which has been exploited to substantially 
develop our understanding of jamming transitions and phase separation phenomena. 
In the framework of zero-range processes, these phenomena may be understood in 
terms of condensation — the property that a single site may acquire a finite fraction 
of the total mass in the system. Moreover, the steady state assumes a relatively simple 
form, namely, a factorised form, in which the steady state is expressed as a product 
of factors, one factor for each site of the system. 

For models with a single conservation law and mass transfer rules which depend 
only on the details of the departure site, the conditions on the hop rates under which 
the steady state factorises have now been established rather completely . However 
systems with many conservation laws are also of wide interest J2| ■ In the soluble cases 
identified so far, the inclusion of additional conservation laws has been shown to enrich 
the phase behaviour dramatically and introduces novel mechanisms of condensation 
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Practical applications of these mass transfer models with multiple conservation 
laws include shaken polydisperse granular systems ^1] and directed networks ^3 E| • 
It is the purpose of this article to derive conditions on the mass transfer dynamics 
under which the steady state assumes a factorised form for a very wide class of models 
with an arbitrary number Q of conservation laws. In section [5] we define a one- 
dimensional, totally asymmetric model with Q species of conserved mass and present 
the factorised steady state. The proof and the necessary and sufficient conditions 
under which it holds are given in section The model is defined with discrete-time 
parallel update, which contains the random sequential update and continuous time 
limits shown in section^] The generalisation to a Q-species mass transfer model on an 
arbitrary lattice is presented in section [5] where sufficient conditions for factorisation 
of the steady state are stated. A number of models which represent applications of 
the steady states derived here are provided throughout. 

2. Q-species mass transfer model in one dimension 

The model is defined on a periodic chain containing L sites. Q species of masses reside 
on the lattice. Associated with each site i is a mass vector mi, the elements of which 

represent the mass of each species q = 1,...,Q present at i: tjh = (m^ , . . . , m$ ). 

The mass of species q at site i, m[ q \ is a continuous variable. 

The dynamics are defined in discrete time, such that the mass vectors at each 
site are updated in parallel during each time step, and the mass is transferred to the 
nearest neighbour site to the right. The mass of species q transferred from site i to 
site i+1 during the time step is /i,- 9 '. The vector (j,. = (fi\ l \ . . . , /4*^) represents the 
masses of all species which move from i to i+1. 

The definition of the model is completed by defining the dynamics. These 
are specified through the mass transfer probabilities tl)i(fi.\mf) which determine the 

stochastic variables fi\ q \ Hence, the probability of transferring the masses /i. from i 
to i+1 in a time step depends on the mass of each species at the departure site, m i; 
and on the site of departure (through the i subscript on ip) . The dynamics conserve 
to total mass in the system of each species. Further, the interactions are zero-range 
— they depend only on details of the site of departure of the mass. 

2.1. Steady state 

The aim is to establish the conditions under which steady state of the model defined 
above assumes a factorised form. Therefore we look for conditions on the choice of the 
mass transfer probabilities under which the steady state weight F^m^) to observe 
the system in a configuration {mj = m 1 , . . . , m L is given by 

L 

F({rn t }) = l[f l (m i ) , (1) 

i=i 

i.e., one factor for each site of the system. We find a necessary and sufficient condition 
for factorisation is that one can find functions and Wi{rn — /i) such that the mass 
transfer probabilities can be expressed in the form 
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where the * denotes a convolution (see l |14[l ). If the mass transfer probabilities are of 
the form ©, then the steady state factorises with single-site weights ^(raj given by 

fi{m.i) = [v * Wi] (mj . (3) 



3. Proof of the steady state 

In this section, we turn to the proof of the steady state to (JSJ . The proof follows 
a similar argument to that given in ^| ^] for the single species case. We write 
F({TOj}, t) to represent the weight of a configuration at time t. For a single update the 
evolution of F^m^}, t), under the dynamics defined by the mass transfer probabilities, 
is described by a master equation which may be written 

L 



■Mi / d Q (j,Mn.\m i ) 



(4) 



where d^m^ represents an integration over the masses of all the species at site i: 
d Q ™; = Uq=i d™i 9) ; similarly, the notation S Q (x) = ]J^ =1 5(x^). 

In the steady state, F({mJ,t+l) = ^({mj,*) = ^({mj). Making these 
replacements in Q and assuming ^({m^}) is given by a factorised form yields 



n/*K) = n / d<3 ^ / ^M^ifeimi) 

2—1 i L 



(•5) 



x ^ Q (™'i - 2Zti + - M^) fi(Mi) ■ 

The aim of the remainder of this section is to find the conditions under which the 
assumption of factorisation holds, i.e., the conditions under which we can solve (O- 
To proceed, we introduce the Laplace transform 



'Mm) 



where s, = (s\ , . 



, s\ Q >) and s ■ m = J2 q s (g W< Hence 

U.M = II / dQ ^J & Q nMa^ 



(6) 



(7) 



At this stage, it is convenient to introduce the function T , i(p i ,m i — (X.) defined as 

Pii^m - a) = ^(/fjzsi) fi(mi) , (8) 

then, doing the integral over the mass variables for each species, mf \ and re-ordering 
the sum in the exponential leads to 

flM = II [/ d^e-*^-^-^ Vi^mi-ld ■ (9) 

i—1 i 

Now, the term in the square brackets on the rhs can be identified with the function 



XiQu, §4+i) = / dPne-^-Oj-^B* Vi (/i , m i — fj, ) 



(10) 
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L L 

n^(«i)=n^(&'&+i)> en) 
i=i i=i 

in order to establish the conditions for factorisation of the steady state. Note that 
we have only rewritten (JjJJ at this stage — a necessary and sufficient solution of (|ll|l 
represents a necessary and sufficient condition for factorisation. 

By taking the logarithm of (|11|) and differentiating with respect to Si then Sj+i, 
it can be shown ^2 Id that a necessary and sufficient condition to solve (|llfl is 

X i (s i ,i Li+1 ) = a i (s i )p(s i+1 ) , (12) 

where the functions CKi(s) and /3(s) are to be determined. This implies 

MSi)=«i(Si)/?(fii)- (1 3 ) 

Since this is just a product over Laplace transforms, its inverse (i.e., fi(m.i)) is a 
convolution 

Mm) = [v* Wi] (m) = / d Q /i u(/i) Wi(m - fi) , (14) 



as given in J2J. 

The functions cxi(s) and /3(s) are determined by 

= J d Q ue-^ Wi (E) , (15) 

where £ = m — /i, and 

« = / d^e^'^M- (16) 

By substituting these into l|12|l . and using (|10|l . one can read off 

Vii^Wn -/£.) = v(fi.)w l {m 1 - fi), (17) 
then from the definition of given in JHJ, equation J21 immediately follows. 

4. Continuous time limit 

The discrete time update contains as a special case random sequential update which 
in turn contains the continuous time limit. The random sequential limit emerges when 
the mass transport probabilities are proportional to a small time step dt such that at 
most one mass transfer event on the lattice takes place within a time step. Then the 
limit dt — > yields continuous time dynamics where the masses are transferred with 
a rate per unit time. We take the random sequential limit first by redefining 

- & Q {^) +x( t i)dt . (18) 

Upon substitution into this yields 

/i(ffii) = ^fc)* Q M + dt[x * Wi]{md + 0(dt 2 ) , (19) 

for the single-site weights, and from J5J) it yields 

MtAm) = 5Q (p) - W .^.) ^ Q (^)[ X * w i](mi) 

i(/i)wi(m,- — n) 
+ dt V - } + 0(dt 2 ) , (20) 
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for the mass transfer probabilities. In the next step, we take the continuous time limit 
by taking dt — ► 0, hence the single-site weights are 

fi(m) = Wi(m) , (21) 

and the hop rates are read off as 

x(u) Wi(m — u) 

«i A* m, = ~ (~ ~ ■ (22) 
— w l {m) 

This hop rate represents the rate at which masses fjS 1 ' , ■ ■ ■ , fjt^' are simultaneously 
transferred from site i to site i + If the hop rates can be written in the form (11'2[l . 
the steady state factorises with single-site weights given by (21). 

We note that the form (1221) implies a constraint on the choice of hop rates. Viewed 
as a recursion, one can express Wj (m) in terms either Wj (m — fi) or Wi(m — fjf), before a 
second application of the recursion yields an expression in terms of Wi(m — u — //); this 
final expression must be the same regardless of the intermediate step, which demands 
that the hop rates satisfy 

Ui((Am) _ Ui(jj/\m) 

Ui(fi\m - fj/) Ui(jj/\m- /£) 

for all i,m, [i and //. This constraint is not an additional constraint to (|22|l : if the 
hop rates can be written in the form (|22|l . they automatically satisfy However, 
it provides an alternative test for factorisation for a given set of hop rates. Moreover, 
this constraint represents a necessary and sufficient condition for factorisation. 



4-1- Two-species ZRP with disorder 

The ZRP with two species of particles is a discrete-mass model in which single units 
of mass, either of species 1 or of species 2, are transferred to the nearest neighbour 
site to the right EH] . Here, the model is generalised to include disorder in the hop 
rates. In the current notation, hop rates for species 1 and species 2 are, respectively, 

^(i,0K (1 Uf)) = ^f ; ^-^) , (24) 

fi(m\ ,m\ ') 

r n ii (i) (2h z(0,l) fj{mf\mf ] - 1) 

Ui{0,l\my,ml J ) = , (25) 

fi(m\ ',m\ >) 

which imply the constraint 

Ul (1 , 0|m t a) , mf ] ) Ui (0, 1 1 mf ) , mf ] ) 



UiCL.Olm^.mf - 1) Ui(0,l|tn, (1) - l,m 4 (2) ) 



(26) 



The single-site weights, obtained by iterating 1)24125(1 . are 

m (1 > m< 2 > 

/i(m (1) ,m (2) ) - n u i (l,0\x,m ( ' 2) )- 1 ]J ^(0, 1|0, y)' 1 , (27) 



having set fi(0, 0) = 1 and absorbed cc(l,0) and a;(0, 1) into the definitions of the 
hop rates without loss of generality. So, for any given hop rates, Il2(it represents a 
necessary and sufficient condition for the steady state to factorise with f,; (m) given by 
H27J1 . This generalises the results of |19l I2TT] to the disordered case. 
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4-2. A two-species model 

In this model, again wc consider two species of discrete masses. Now, the number of 
particles of a species which hops depends on the number of particles of the other species 
at the departure site. Hence, we consider the following rates: if rn^p > rnp > 0, then 

(2) 

ml particles of species 1 hop with a rate 

f(m\ \m\ ') 

if m\ > mp > 0, then mp particles of species 2 hop with a rate 

»(0^ ll, rf , ,^)=^^ ! f##i (29) 

f(m\ ,m\>) 

if 77ij = rn'p = m > then, with equal probability, either m particles of species 1 
hop with a rate 

u(m,0\m,m) = f,^' m , , (30) 
f{m,m) 

or m particles of species 2 hop with a rate 

/n i ^ f(m,0) ( . 

u({j,m\m,m) = —. r . (31) 

f (m, m) 

If only a single species is present at a site, single particles hop with rates 

*°Koh«^>, *.M-^2. (-) 

A simple choice for the rates is u(n, 0|m, n) = (1 + n/m) b for m > n > and 
u(0,m\m,n) = (1 + m/n) b for n > m > 0; the single-site weights for these rates 
are given by f(m, n) = (m + n)~ b . Though the model presented here is new, single- 
site weights of this form have been considered in |16| where it was shown that a 
condensation transition occurs above a critical particle density for b > 3. Ordinarily, 
the two species condense at completely independent sites; in the present model 
however, an exclusion interaction exists between the condensates: by considering the 
hop rates in the limit of large n and m, it is clear that if both species condense at 
the same time, they do so at randomly located sites subject to the condition that the 
condensates of each species do not occupy the same site. 



5. Generalisation to hypercubic and arbitrary lattices 

The models considered above can be generalised to Q-species mass-transfer models 
on arbitrary lattices. The derivation of the conditions for factorisation can also be 
generalised along the lines of 11 . Though this derivation is straightforward, its 
simplicity is rather obscured under complicated notation; therefore we present in this 
section the results without providing details. We also provide concrete examples of 
models with factorised steady states in order to illustrate the utility of the results. 
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5.1. Generalised model 

Again, we begin with discrete time dynamics whereby each site is updated in parallel. 
The new ingredient is that Q fixed sets of directed links now connect certain pairs of 
sites — each set is associated with a particular species and mass of any species can 
only be transferred along associated directed links. 

To define the update rules, we define a mass transfer matrix for each species q 
with elements (x^ . These elements represent the mass of species q which is transferred 
from site i to site j during a time step. The vector . = (fJy \ ■ ■ ■ , A 1 ^) represents 
the masses of all species which move from i to j] if there is no directed link for species 
q pointing from i to j, then = identically. Also, fj,^ is in general nonzero and 
represents the mass of species q at site i that does not hop during the time-step. 

The elements \xf^ which are not set to zero for all times (i.e., those which are 
accompanied by a directed link pointing from i to j) are determined by the mass 
transfer probabilities ^({/i J-lmJ. The set {a*^} contains all sites j to which i is 
connected by a directed link for at least one species and includes site i itself. Thus it 
specifies the mass transferred for each and every species which can hop from i to j. 

As before, the dynamics conserve the total mass of each species in the system. 
The sum over the elements in a row of the mass transfer matrix for species q is equal 
to the mass of species q at the site before the update: J2j Inf — TO ! 9 ' '■ Similarly, the 
sum over column elements represents the mass of species q at i after the update: 



(?) _ '(<?) 



rn ,■ 



5.2. Conditions under which the steady state factorises 

Using the multi-species procedure presented in section [3] and following the technique 
outlined in JT] for one species of mass on an arbitrary lattice, one can establish 
conditions on the generalised mass transfer probabilities under which the steady state 
assumes the factorised form (Q. 

We find a sufficient condition for factorisation (but no longer necessary in general) 
is that one can find functions Vij(p i .) such that the mass transfer probabilities can be 
expressed in the form 

11, 

Mii*J\sk) = t- 1 — fr— . ( 33 ) 



-y 



n*i v ji (m) 



subject to the constraint J2j f^if = m ^ f° r au species q. Here, the product in the 
numerator is over sites j which have a directed link coming from site i for at least 
one species; the product in the denominator is over all sites which have at least one 
directed link pointing into i and the * denotes a convolution form. If the mass transfer 
probabilities are of the form JSS}, then the steady state factorises with single-site 
weights fi{rrij) given by 



fiinii) 



n 



*3 



(mi) . (34) 



which again is a multiple convolution over those sites which have a directed link into 
i for at least one species. 
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A second condition for factorisation is that the functions (/•*••) satisfy 



Yl v a 


(Mi) = 




(mi) , 


(35) 


*j 




*j 







which emerges as a constraint on the geometry of the underlying lattice. This 
constraint will typically limit the choices allowed for the dynamics on inhomogeneous 
lattices. On homogeneous lattices, there exist two particular classes of dynamics 
for which the constraint is always satisfied: (i) detailed balance — in cases where 
the dynamics are symmetric, i.e., Vij(fX..) — Vji(ji..) the constraint (|35fl is always 
satisfied; (ii) pair-wise detailed balance — the constraint H35|l is always satisfied in 
cases where, for each species, there exists a site i with a link into j for every link from 
j to another site k, such that Vij(fj,._.) — Vjk(^- k )- A simple example is asymmetric 
nearest neighbour mass transport in one dimension, where transfer to i from i — 1 
can be paired with transfer from i to i + 1. There are many more ways to satisfy the 
constraint, some of which are discussed below, but dynamics which fall into either of 
the above two categories are sufficiently prevalent in mass transport models that we 
wish to emphasise that the constraint is automatically fulfilled in these two cases. 

The general conditions under which (|33ll represents a necessary condition seem 
difficult to formulate [TT]: though necessary for totally asymmetric mass transport in 
one dimension, as outlined above, it is not necessary for symmetric mass transport. 
It is also necessary for mass transfer on the complete graph. 



5.3. Continuous time limit 

The random sequential then continuous time limits are obtained as before, by 
redefining 

% (Hij ) = sQ (Hij ) + x v (My)di ■ (36) 
In the limit dt — > 0, the single-site weights assume the form 

fi(md = v u (Mi) > ( 37 ) 
and hop rates are obtained as 

Xij (/i) va ( m i — [i) 

UijWmj = — = — t — v — — • ( 38 ) 

- vu(mi) 

This hop rate represents the rate at which masses ■ ■ ■ , are simultaneously 
transferred from site i to site j. Finally, in the continuous time limit, the constraint 
acquires the form 

EN * = * w «K™i) • (39) 

By taking the Laplace transform, cancelling common factors then inverting Laplace 
transform back, the constraint can be written 

Therefore, if the hop rates can be written in the form and provided the constraint 
l|4U|l is satisfied, the steady state factorises with single-site weights given by l|37fl . 
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5.4- Generalised Q -species mass-transfer model in one dimension 

In the case where each particle species is transferred between nearest neighbour sites 
on the same one-dimensional lattice, the constraint 14U|) implies 

£'M+i(m) +Xi,i-i(fj) = Xi+i,»(/£) +Xi-i ti (fi) , (41) 

and the general solution has been obtained in [23 ; this solution must be used in order 
to allow disorder in the asymmetry. Here, we take the pair-wise balance solution 

Xi,i+x{(j) =p(n) , (42) 
Xi+i,i(f±) = q([i) , (43) 
for all i. Hence the steady state factorises if the hop rates can be written in the form 

uu+iitAmJ = — — , (44) 

, q(j±) fiiMi-tj) ,._v 

«*.<-itelm*) = Jt^ • ( 45 ) 

for all i. Note that some disorder is still incorporated through the site dependence in 
Vu (m) : the disorder is in the departure rates rather than the asymmetry. 

The two-species, totally asymmetric limit q(fJ>) = with transfer of single units 
of discrete mass recovers the limit discussed in section IQ1 

5.5. Generalised Q-species mass transfer on a hypercubic lattice 

Here, we consider homogeneous, symmetric dynamics on a hypercubic lattice: masses 
/i = (/Z 1 ), . . . , M ) are transferred to a randomly selected nearest neighbour site with 
a rate u(p\m). 

In this case, the dynamics satisfy detailed balance so the constraint H4Q(I is easily 
seen to be satisfied for Xij(fi) independent of i and j. The sufficient condition for 
factorisation then is that the hop rates satisfy the constraint 

u(tAmi) _ uifj/lrrii) 

u{tAmi~li) u^'lmj-^) 5 

in which case the single-site weights are determined from 

f(m i — fi) 

UWXlli) = —77 — 7^ , (47) 

where the constant factors x(fi) have been absorbed into u. Thus, since the rates 
are homogeneous, the single-site weights are independent of i. This generalises the 
results of [22] to an arbitrary number of species of particles. Generalisations to include 
disorder or to partial asymmetry along the lines of the section[^31 are straightforward. 

5.6. Two-species ZRP on an arbitrary lattice 

In this section, we consider the two-species ZRP again, but this time on an arbitrary 
lattice. We consider the case where both species hop on the same underlying lattice, 
with hop rates from site i to site j given by 

u y (l > 0|mj 1) ,mf ) ) =y ij ^(m^mf ) , (48) 
Uij (0, l\m?\r4 2) ) = % af\mf\mf ] ) , (49) 
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for species 1 and 2, respectively. Here, af^ (m^~\ rn^ 1 ) represents the total departure 
rate of particles of species q from site i, for q = 1,2; represents the probability that 
a particle is transferred from site i to site j. The same probabilities are used for both 
species because they are both moving on the same underlying lattice (though this is 
straightforward to generalise). 
The choices 

Xij (m) = VijPi , (50) 
for /i = (1, 0) or (0, 1), and 

1 1 1 ,m ' J ~ / 4 (m(D,m(2)) ' 
are of the form (|38|l from which the single-site weights are obtained as 

m (1) m (2) 

/,(!»« =p^ 1 ^ II «! 1) (- I - (2) )" 1 II ^(O,*)" 1 ■ (53) 

a:— 1 y— 1 

From the constraint H40|) . the pj's are determined by 

Pi = y i ip i ' ( 54 ) 

which is just the solution for the steady state of a single random walker moving on 
the lattice with hopping probabilities defined by the ytj (c.f. I0TJ). 

Again the two equations (|51l I52|) imply a constraint on the possible choices of 
hop rates for which the steady state still factorises: 

(1)/ (1) (2K (2)/ (1) (2h 



(2)/ (1) (2) ,l (2)/ (1) -, (2)n ' 

If this equation is satisfied, a model with rates (|48l I49fl factorises with single-site 
weights given by with the p^s determined by — thus it is the p^s which 
contain the information about the structure of the underlying lattice. 



6. Conclusion 



We have found conditions for factorised steady states applicable to a very wide class 
of models. The results have been exploited to write down the exact steady states 
for several models which generalise previously solved cases. The interplay between 
conservation laws, geometry, disorder and the nature of interactions can all be explored 
within the framework of the steady states derived here and merit further investigation. 

As intimated in the introduction, condensation is of particular interest in mass 
transfer models. For homogeneous two-species systems, the general conditions for 
condensation are well understood at the level of the steady state [23] - However the 
models considered here are substantially more general, therefore it is of interest to 
establish whether new condensed phases emerge. Even for the two species model 
considered in section FOl which has a steady state in the class considered in [55], the 
exclusion interaction between condensates seems to be a property of the dynamics 
which has not been observed in previously studied models. 
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Moreover, the coarsening dynamics in two-species models have been studied using 
heuristic, quasi-static scaling arguments |2H, and are far richer than their single- 
species counterparts. It is reasonable to suppose that the wider class of models solved 
here will further enrich the phenomenology of the coarsening dynamics. 

Finally, the question of the structure of steady states beyond the factorised form 
remains open |25| but must be addressed in order to understand behaviour observed 
in non-factorisable models |2fi| . 
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